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Abstract 

We analyse the production and propagation of pions and kaons in heavy-ion re- 
actions from 0.8 - 1.8 A-GeV within a coupled channel transport approach including 
the kaon production channels BB -> K + YN, irB -> K + Y, BB -» NNKK, irB -» 
NKK, K + B — > K + B and tttt — * KK. Assuming the hyperon selfenergy to be 2/3 
of the nucleon selfenergy we find that all inclusive experimental K + spectra at SIS 
energies can be reproduced reasonably well without any selfenergies for the kaons 
although a slightly repulsive kaon potential cannot be excluded by the present data 
on kaon spectra and flow. 
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1 Introduction 



The study of hadron properties in the dense nuclear medium via relativistic nucleus- 
nucleus collisions is a major aim of high energy heavy-ion physics. Especially the pro- 
duction of particles at 'subtreshold' energies is expected to provide some valuable insight 
into the properties of hadrons at high baryon density and temperature |jj . In this respect 
selfenergy effects in the production of particles have been found previously for antiprotons 
by a number of groups @-[[| though the actual magnitude of the attractive p-potential in 
the nuclear medium is still a matter of debate. Also antikaons according to Refs. |6|-||15|| 
should feel strong attractive forces in the medium whereas the kaon potential is expected 
to be slightly repulsive at finite baryon density. A first exploratory study with respect 
to antikaon selfenergies by Li, Ko and Fang in Ref. [pl| indicated that sizeable attractive 
K~ potentials are needed to explain the experimental spectra from [|17}] for Ni + Ni at 
1.85 A-GeV . Indeed, their findings could be substantiated recently in Ref. |T!| in a system- 
atic analysis of antikaon production in nucleus-nucleus collisions. So far the production 
of K + mesons has been primarily addressed in order to learn about the incompressibility 



of the nuclear equation of state p9| , pOf . On the other hand, this process could also give 
information on the predicted repulsive kaon potential. This should also be seen in reduced 
kaon production yields since at 1 - 2 A-GeV already considerable baryon densities (up to 
3 po ~ 0.5 fm -3 ) can be reached in the compression phase According to the study in 
Ref. 15j the repulsive kaon potential might be even about 100 MeV at 2po which should 
lead to a sizeable reduction of the cross section. 

The most serious problem related to K + production at 'subthreshold' energies are the 
baryon-baryon elementary production cross sections close to threshold where no experi- 
mental data had been available so far and rough extrapolations have been used PTl 
Only very recently a new experimental data point has been obtained 2 MeV above the 
kaon production threshold in pp collisions p3| , which can be used for more accurate 
parametrizations of the elementary production cross section and a systematic reanalysis 
of the experimental data taken up to now. 

Our paper is organized as follows: In Section 2 we briefly describe the transport 
approach employed as well as the new parametrizations for the elementary strangeness 
production channels. Section 3 contains a detailed comparison of our calculations with 
the available data at SIS energies while Section 4 concludes our study with a summary. 
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2 Model — Ingredients 



In this paper we perform our analysis along the line of the HSD[] approach |24|] which is 
based on a coupled set of covariant transport equations for the phase-space distributions 
fh(x,p) of hadron h |24], |25| 



i.e. 



{(n, - n v %u;; - M* h d;u s h ) an + (n u d*u» + M* h d*u s h ) &>;} f h (x, P ) 

= J2 J d2d3dA...[G^G] 12 ^ 34 ..5\n + U 2 -n 3 -U A ...) 

h2hzh,4,... 

x { fh 3 (x, Vi)fh< (x,p 4 ) f h (x, p) f h2 (x, p 2 ) 

- fh(x,p)fh 2 (x,P2)fh 3 {x,P3)f h4 (x,p 4 )} . . • • (1) 

In Eq. (H) U^(x,p) and Uj^(x,p) denote the real part of the scalar and vector hadron 
selfenergies, respectively, while 12^34... (^(II + Ii 2 — n 3 — n 4 . . .) is the 'transition 

rate' for the process 1 + 2— > 3 + 4 + . . . which is taken to be on-shell in the semiclassical 
limit adopted. The hadron quasi-particle properties in (P are defined via the mass-shell 
constraint p5| , 

5(U,W-M* h 2 ) , (2) 
with effective masses and momenta (in local Thomas-Fermi approximation) given by 

M*(x,p) = M h + U*(x,p) 

W(x,p) = pP-U£(x,p) , (3) 

while the phase-space factors 

f h (x,p) = l±f h (x,p) (4) 

are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depending 
on the type of hadron in the final/initial channel. The dots in Eq. ([!]) stand for fur- 
ther contributions to the collision term with more than two hadrons in the final/initial 
channels. The transport approach ([I]) is fully specified by U^(x,p) and U£(x,p) (/i = 
0, 1,2,3), which determine the mean-field propagation of the hadrons, and by the tran- 
sition rates G^GS 4 (. . .) in the collision term, that describes the scattering and hadron 
production/absorption rates. 

The scalar and vector mean fields and U£ for baryons are taken from Ref. 



the hyperon selfenergies are assumed to be 2/3 of the nucleon selfenergies as in |[24| . In 
the present approach we propagate explicitly pions, 77's, p's, a/s and <E>'s as free particles 
whereas kaons and antikaons are propagated with effective potentials. 



1 Hadron String Dynamics 
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As in case of antiprotons there are several models for the kaon and antikaon selfenergies 
P HI15H ) which all differ in the actual magnitude of the selfenergies, but agree on the 



relative signs for kaons and antikaons. Thus in line with the kaon-nucleon scattering 
amplitude the K + potential should be slightly repulsive at finite baryon density whereas 
the antikaon should see an attractive potential in the nuclear medium. Without going 
into a detailed discussion of the various models we adopt the more practical point of view, 
that the actual K + selfenergies are unknown and as a guide for our analysis use a linear 
extrapolation of the form (cf. |J 



m K(pB) = m° K ( 1 + a— ) , (5) 

with a describing the strength of the kaon potential at finite baryon density p#. For our 
following analysis we will restrict to a = and a = 0.06 to model a slightly repulsive kaon 
potential. We note that the kaon potential determined from the kaon-nucleon scattering 



length clkn within the impulse approximation [26| leads to a ~ 0.06 when using the 
isospin averaged scattering length cl^n ~ —0.255 fm from Ref. [E7 . 

First, the individual production channels of the kaons have to be specified. Here, we 
express the cross sections as a function of the scaled invariant energy squared sq/s, since 
the change of the quasi-particle mass then can be incorporated in the threshold energy 
y/so for the particular channel. This recipe might be still a matter of debate; our findings 
in Refs. [^8], [29|, |30], |3l|] indicate, that the production is essentially dominated by phase 
space close to threshold and thus a scaling in sq/s should be a good approximation. 

The isospin averaged production cross section of a K + A and K + Yi pair in a nucleon- 
nucleon collision is related to the measured isospin channels as: 

cnn^k+an = 2 a pp-* K+A P 

(6) 

3 

E0 P ) • (7) 



Following |29| the reaction cross section can be approximated by 
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with ^/soi = m^—m^+m^ and ^fsm = m-^+m^+m^. According to isospin relations the 
NA and A A production channels get additional factors of 3/4 and 1/2, respectively. This 
scaling of the AiV and A A production channels from Ref. ||22|| , however, is questionable 
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and further microscopic studies along the lines of Refs. |32|, |33| will be necessary to 
determine these reaction cross sections more accurately. 

The elementary cross sections for the pion induced channels ttN — > K + Y have been 
computed by Tsushima et al. in Ref. [34], which we adopt for our present study. K + 
elastic scattering with nucleons also has an impact on the final kaon spectra; the elas- 
tic cross section employed is displayed in Fig. 5 of Ref. HIBfl . In addition, we include 
the KK production channels in baryon-baryon {BB), ttB and tttt collisions within the 
parametrizations from |1| which, however, do not contribute significantly to the inclusive 
K + yield. 

The calculation of 'subthreshold' particle production is described in detail in Refs. |], 
35p and has to be treated perturbatively in the energy regime of interest here due to the 
small cross sections involved. Since we work within the parallel ensemble algorithm, each 
parallel run of the transport calculation can be considered approximately as an individual 
reaction event, where binary reactions in the entrance channel at given invariant energy 
a/s lead to final states with 2 (e.g. K + Y in ttB channels) or 3 (e.g. for K + YN channels 
in BB collisions) particles with a relative weight Wi for each event i which is defined by 
the ratio of the production cross section to the total hadron-hadron cross section]]. The 
perturbative treatment now implies that in case of strangeness production channels the 
initial hadrons are not modified in the respective final channel. On the other hand, each 
strange hadron is represented by a testparticle with weight Wi and propagated according 
to the Hamilton equations of motion. Elastic and inelastic reactions with pions, r/'s or 
nonstrange baryons are computed in the standard way however, only the dynamical 
feedback of the strange hadrons to the nonstrange mesons and baryons is neglected. The 
final cross section is obtained by multiplying each testparticle with its weight In this 
way one achieves a time-saving simulation of the strangeness production, propagation and 
reabsorption during the heavy-ion collision. 



3 Results 

Since A- and pion-induced channels play a major role for K + production, we start our 
analysis with a comparison of our calculations for pion production at SIS energies with 
the available n + data for Ni + Ni at 1.0 and 1.8 A-GeV and Au + Au at 1.0 A-GeV 
in Fig. [I]. The experimental spectra at 6^ = 44° ± 4° []36|] are described reasonably 



well in the whole kinematical range not only for Ni + Ni, but also for Au + Au. We 
slightly underestimate the pion spectrum for Ni + Ni at 1.8 A-GeV at low momenta which 



2 The actual final states are chosen by Monte Carlo sampling according to the 2, or 3-body phase 
space. 
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might reflect limitations of our configuration space showing up at higher bombarding 
energy. For future comparison we also include our results for the inclusive tt + spectra 
from Au + Au at 1.5 A-GeV (upper histogram). For a more detailed analysis of pion 
production at SIS energies we refer the reader to Ref. |37 |, which works in a very similar 
theoretical framework. For our present purpose we only conclude that the pion and baryon 
dynamics including the A and higher resonance dynamics in our transport calculations 
are sufficiently well under control. 

The Lorentz invariant K + spectra for Ni + Ni at 0.8, 1.0 and 1.8 A-GeV from the 
same calculations as above are shown in Fig. ^|in comparison to the data from | 3"8fl . Here 
the full lines reflect calculations including only bare K + masses (a = 0) while the dashed 
lines correspond to calculations with a = 0.06 in Eq. ([5]), which leads to about a 30 MeV 



increase of the kaon mass at po- Note that the repulsive kaon potential from Ref. |15| is 
best fitted for a = 0.1, which would imply a further reduction of the K + cross section. 
When fitting an exponential to the Lorentz invariant spectrum ~ exp(— E/T Q ), where E 
is the kinetic energy of the kaon in the cms, we obtain slope parameters T ~ 60 MeV, 
75 MeV and 102 MeV at the bombarding energy of 0.8, 1.0 and 1.8 A-GeV , respectively. 

The general tendency seen at all bombarding energies is that our calculations with a 
bare kaon mass provide a better description of the experimental data for Ni + Ni than 
those with an enhanced kaon mass. Since this general trend might be accidental we 
also compare our calculations for the heavier systems, i.e. Bi + Pb at 0.8 A-GeV and 
Au + Au at 1.0 A-GeV , with the respective experimental data [11]] in Fig. |3]. The 

full dots represent the earlier data for Au + Au at 1.0 A-GeV [EOI while the open circles 



result from a new measurement |4lJ of the system at the same bombarding energy. Both 
calculations (solid line: a = 0; dashed line: a = 0.06) are compatible with the data due to 
the experimental uncertainties. As in case of antikaons [18[] the selfenergy effects are most 
pronounced at low center-of-mass momenta of the kaons where, unfortunately, no data 
have been taken so far. For future comparison we also include our results for Au + Au 
at 1.5 A-GeV . 

With our new elementary cross sections the relative weights of the K + production 
channels change considerably compared to our earlier calculations |32[] as can be seen 
from Fig. [| for Au + Au at 1 A-GeV . Whereas the NN production channels (dashed 
line) are almost negligible as before [[|2|] the dominant yield now stems from irN reactions 
(dot-dot-dashed line) which surpass the iVA channel (dot-dashed line). This change 



comes about because our old calculations |42] were performed in the so-called frozen-A 
approximation in which the A's are not allowed to decay until the collision has ended. 
The main message of the old studies, that the secondary reactions are very important 
for a description of the kaon yield, however, survives. The calculation in Fig. El has been 
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performed for a bare K + mass; we note that the relative channel decomposition does not 
change very much when performing a calculation with a slightly repulsive kaon potential. 

In Fig. |5| we present the ratio Nx+/N n + versus the number of participating nucleons 
Apart - which is a measure of the impact parameter b - from our calculations for Ni + Ni 
at 1.0 and 1.8 A-GeV as well as for Au + Au at 1.0 and 1.5 A-GeV for a = (solid lines) 
and a = 0.06 (dashed lines). The experimental data have been taken from Ref. |[4U|| . The 
'theoretical' error bars indicate the calculational uncertainty due to the finite particle 
statistics, since only mesons for 9^ = 44° ± 4° have been taken into account. The ratio 
Nx+/N n + increases with A part in all cases; this increase is more pronounced at the lower 
bombarding energy of 1.0 A-GeV and for the heavier system Au + Au. In our calculation 
this increase in the relative K + yield is due to a strong rise in N K +/A part ; N n /A part is flat 
as a function of A part . Because the kaons are mainly due to the secondary collisions, they 
are predominantly produced at high baryon density 1.5po < Pb < 2.5po while the pion 
production occurs at lower density, too. Since the volume for the high baryon density 
increases sizeably when going from peripheral to more central collisions, the K + /ir + ratio 
has to increase accordingly. The relative contribution of the various production channels 
to the total K + yield are shown in Fig. ^ as a function of A part for Ni + Ni at 1.0 and 
1.8 A-GeV . It can clearly be seen that the relative importance of the primary AW channel 
increases with energy, but that of the secondary ttN reaction decreases. This explains the 
flatter behavior seen in Fig. [| for the higher energies. 

In Fig. |7| we show the angular distribution of the produced kaons for Au + Au at 
1 A-GeV which exhibits a forward-backward peaking. By comparing our calculations with 
and without angular dependence in the elementary kaon production cross section we have 
verified that this angular anisotropy is not due to the elementary production processes. 



Instead, it reflects the forward-backward peaking of the pion angular distribution |37 



The figure shows that the rescattering of the kaons enhances the anisotropy somewhat, 
but is not solely responsible for it. 

A further quantity of interest - in the context of our present analysis - is the kaon 
flow in the reaction plane, which should show some sensitivity to the kaon potential in 



the nuclear medium as put forward by Li, Ko and Brown p3| , |4"4" , [45| , f46|| and is presently 
investigated by the FOPI Collaboration [|7 , l"8fl . Here due to elastic scattering with 
nucleons the kaons partly flow in the direction of the nucleons thus showing a positive 
flow in case of no mean- field potentials P3|| . With increasing repulsive kaon potential the 
positive flow will turn to zero and then become negative; experimental data on kaon flow 
thus are expected to discriminate further on the potentials seen in the medium. In order to 
investigate this question we have performed detailed (and high statistics) calculations for 
K + production in Ni + Ni reactions at 1.93 A-GeV as measured by the FOPI Collaboration 
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48 1 - In order to compare with their data we have included a transverse momentum cut 



Pt > 0.5 rriK, where is the kaon mass, and gated on central collisions with impact 



parameter b < 4 fm as in Ref. []4"5 ]. The results of our calculations are displayed in Fig. 
|] in terms of < p x > / itlk versus the normalized rapidity 

y° = -^, (11) 

Uproj 

where y cm and y pro j are the kaon and projectile rapidity in the cms, respectively. Our 
calculations without any kaon potential (a = 0, full line) indeed show a positive flow 
as expected, which appears still to be compatible with the data within the error bars 
and is practically identical to the respective calculations from Ref. [45]. On the other 
hand, increasing the kaon potential (a = 0.06, dashed line) the kaon flow becomes slightly 
negative or comparable to zero, quantitatively in line with the calculations from Ref. E5 



and in somewhat better agreement with the data. Thus in case of the flow observable a 
slightly repulsive kaon potential (~ 30 MeV at po) is more in line with the present data 
of the FOPI Collaboration M . 



4 Summary 

In this work we have presented a detailed study of pion and kaon production in nucleus- 
nucleus collisions for medium and heavy systems from 0.8 to 1.8 A-GeV within the coupled 
channel BUU approach, where the kaons are produced perturbatively, however, propa- 
gated explicitly with their final state interactions. An important ingredient of our reanal- 
ysis of the kaon cross sections are the novel elementary production cross sections from 



Refs. [28 , [29| , |3l|| and from Ref. |34| for the pion induced channels. We note, however, that 
the AN and AA production channels are not that well determined and require further 
theoretical efforts. 

Our analysis shows that the tt + spectra are reasonably well described in this energy 
regime without introducing any medium modifications for these mesons (cf. also Ref. |37|). 
This is understood in terms of the strong reabsorption of pions that essentially leads to 
the surface-emission. For the much more penetrating K + mesons our results also show no 
final convincing indications of in-medium changes. The new kaon spectra for the heavy 
system Au + Au at 1 A-GeV as well as the kaon flow data for Ni + Ni at 1.93 A-GeV can 
- within their errorbars - be described by assuming free kaon properties although there 
seems to be a tendency to favor a small repulsive kaon potential of about +30 MeV at 
normal nuclear matter density, which would be in line with the kaon potential as extracted 
from the kaon-nucleon scattering length in the impulse approximation PS |. On the other 



hand, we find that more repulsive kaon potentials as predicted by some Lagrangian models 
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appear not to be compatible with the available data and in particular with the K + /n + 
ratios. An upcoming data analysis in the form of Fig. [5] and high statistics data on kaon 
flow in the form of Fig. [8] should shed some further light on this issue. However, even 
then theoretical uncertainties connected with our insufficient knowledge of A-induced 
cross sections will remain. 

The authors acknowledge valuable and inspiring discussions throughout this work with 
N. Herrmann, H. Oeschler, P. Senger, A. Sibirtsev and Gy. Wolf. 
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Figure 1: Inclusive 7r + spectra from Ni + Ni collisions at 1.0 and 1.8 A-GeV and Au + Au 
at 1.0 A-GeV in comparison to the experimental data from Ref. |3(| at 6i a b = 44° ± 4°. 
The upper solid histogram displays our calculations for Au + Au at 1.5 A-GeV . 
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Figure 2: Inclusive K + spectra from Ni + Ni collisions at 0.8, 1.0 and 1.8 A-GeV in 
comparison to the experimental data from Ref. at 9i a b = 44° ± 4°. The solid lines 
represent calculations with bare K + masses, while the dashed lines result for a = 0.06 in 
Eq. (|). 
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Figure 3: Inclusive K + spectra from Bi + Pb collisions at 0.8 A-GeV and Au + Au at 
1.0 A-GeV in comparison to the experimental data from Ref. |M [4(1 Ell at 9^ = 44° ±4°. 



The solid lines represent calculations with bare K + masses, while the dashed lines result 
for a = 0.06 in Eq. ([5]). The calculations for Au + Au at 1.5 A-GeV are included for 
future comparison. 
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Figure 4: Inclusive K + spectra from Au + Au collisions at 1.0 A-GeV (solid line) in 
comparison to the experimental data from Ref. [|40"| , 41] at 6^ = 44° ± 4° for the bare 
kaon mass scenario. The dashed line represents the contribution from NN collisions 
while the dot-dashed and dot-dot-dashed line show the contributions from AN and nN 
collisions, respectively. 
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Figure 5: The K + /tt+ ratio for Ni + Ni at 1.0 and 1.8 A-GeV and Au + Au at 1.0 
and 1.5 A-GeV as a function of the number of participating nucleons A part as defined in 
Ref. P0|. The solid lines represent calculations with bare K + masses, while the dashed 



lines result for a = 0.06 in Eq. (Rl). The experimental data have been taken from Ref. ED 
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Figure 6: The relative contribution of the various production channels to the total K + 
yield as a function of A part for Ni + Ni at 1.0 and 1.8 A-GeV . The solid line represents 
the contribution from txN collisions while the dashed and dot-dashed lines show the 
contributions from AN and iViV collisions, respectively. 
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Figure 7: Angular distribution of kaons in the nucleus-nucleus cms for Au + Au collisions 
at 1 A- GeV including all impact parameters. The dashed line represents a calculation 
without K + N rescattering while the solid line shows the results when including kaon 
rescattering. 
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Figure 8: Kaon flow in the reaction plane (< p x > /nix) as a function of the normalized 
rapidity y° = y C m/y P roj for Ni + Ni at 1.93 A GeV. We have gated on central collisions {b < 
4 fm) and applied a transverse momentum cut for px > 0.5 nix as for the experimental 
data of the FOPI Collaboration |4l| (full dots). The open dots are obtained by reflection 
at y° = 0. The solid line and dashed line display the results of our transport calculations 
without (a = 0) and with a slightly repulsive kaon potential (a = 0.06), respectively. 
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